arXiv:1503.03434v3 [physics.gen-ph] 27Jul2015 


Bicriteria Optimization in Designing Magnetic 
Composites 


Krzysztof Z. Sokalski 

Institute of Computer Science, Czestochowa University of Technology, Al. Armii Krajowej 

17, 4&-200 Czestochowa, Poland 

Barbara Slusarek, Bartosz Jankowski, Marek Przybylski 
Tele and Radio Research Institute, Ratuszowa 11 Street, 03-450 Warszawa, Poland 


Abstract 

A novel algorithm for designing values of technological parameters for produc¬ 
tion of soft magnetic composites (SMC) has been created. These parameters 
are hardening temperature T and compaction pressure p. They enable us to 
optimize power losses and induction. The advantage of the presented algorithm 
consists in bicriteria optimization. The crucial role played by the presented al¬ 
gorithm is the scaling and pseudo-state equation. On this basis mathematical 
models of power losses and induction have been created. The model parameters 
have been calculated on the basis of the power loss characteristics and hystere¬ 
sis loops. The created optimization system has been applied to specimens of 
Somaloy 500. The obtained output consists of a finite set of feasible solutions. 
In order to select a unique solution an additional criterion has been formu¬ 
lated. Keywords: soft magnetic composites; power losses; induction; bicriteria 
optimization 
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1. Introduction 


Soft magnetic composites (SMC) have physical properties which are used 
for adapting these materials to specific applications |l|,|2j|. Very often the func¬ 
tionality of these materials depends on more than one feature. This leads to 
multi-criteria optimization problems, which has not been applied yet in the de¬ 
sign of SMC. However, there are papers which treat more than one physical 
property of SMC but these are not considered as target functions in an opti¬ 
mization procedure ^|,Q{. Recently an algorithm for designing values of the 
hardening temperature and the compaction pressure in the production process 
of soft magnetic composites (SMC) has been derived by using the concept of the 


pseudo-state equation Q. In equilibrium thermodynamics the equation of state 
relates thermodynamic parameters. For instance, in the case of gas-liquid sys¬ 
tem they are the temperature, pressure and volume of the considered material. 
By analogy with the equation of state we consider a phenomenological relation 
between the technological parameters and physical properties of the material. 
Such an approach for SMC is possible thanks to the topology of the completed 
set of scaled power loss characteristics. The most important features of this 


5,Q: 


topology are the following, 

Set of characteristics consists of one variable smooth function 


(B m )P 


= F 


B a 


(1) 


fl: 


where Ptot is density of power loss, B m is peak of induction, / is frequency of 
electromagnetic field wave, F(-) is a function of the following form 

F (Jz) = + f / B m ' F* + f / B m ' ( P 3 + f /• T 4 ))) , (2) 

where Tj, a and 0 have to be determined from experimental data. The form 0 
has been derived from the assumption about the power losses as a homogeneous 
function in a general sense. Each characteristic is determined by the values 
of Tj coefficients and a as well as 0 exponents. These are functions of the 
technological parameters T and p 


Tj = Tj(T,p), a = a(T,p), 0 = 0(T,p), 


( 3 ) 
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where T and p are hardening temperature and compaction pressure, respectively. 
m reveals that characteristics ©, © of samples composed at different T,p 
conditions possess different dimensions, whereas all disentangled characteristics 

Ptot = B&if/BZ • (Id + //B“ • (r 2 + //B“ ■ (r 3 + //B“ ■ r 4 ))), (4) 

possess a common physical dimension. Why do we use the implicit form ©,©? 
Note that the right-hand side of this equation depends only on one effective vari¬ 
able -g^-. Therefore, calculations performed with ©,© are represented by one 
curve for all values of / and B m , whereas results of calculations performed with 
0 are split into many curves. For instance, if one needs Ptot as a function of / 
then the number of generated characteristics is equal to the required number of 
different values of B m . Because of the different dimensions of the different char¬ 
acteristics they do not cross each other except at origin point -g4- = 0, for which 
dimension is not very important. According to the Egenhofer theorem [fj] the re¬ 
lations between characteristics are invariant with respect to scaling, translation 
and rotation. Just the conservation of the relations with respect to the scaling 
enables us to use the implicit form of characteristics. All the characteristics 
are monotonic increasing functions of -Jzr According to m the power loss char¬ 
acteristics are parameterized by pressure and temperature. This dependence 
enables us to introduce a measure of distance in the space of characteristics. 
Let (pi,Ti) and (p 2 • Ti ) be labels of the characteristics of the two composites 
which have been composed under conditions corresponding to these pressures 
and temperatures, respectively. Then the distance between these characteristics 
has the following general form: 

p(p 1 ,T 1 ,p 2 ,T 2 ) = R(\p 1 -p 2 \,\T 1 -T 2 \), (5) 

where R(-, •) satisfies axioms of the distance function. Therefore, the set of all 
characteristics constitutes metric space.We have shown in Q that this space 
consists of two subspaces. Therefore, by introducing the distance measure in 
the space of all characteristics we make each of these subspaces compact. 

Each compact set corresponds to physical phase which is defined by char- 
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Figure 1: The five loss characteristics for T < 550°C corresponding to the low losses phase 
and the one characteristic for T = 600° C corresponding to the high losses phase. 


acteristic values of the physical parameters. For instance, in [5[ we considered 
the low and high losses phases of SOMALOY 500. All the properties mentioned 
above are presented in Fig. [Tj however, the compactness of characteristics’ sub¬ 
sets is ensured by the existence of ©. These properties have enabled us to 
introduce a measure of power loss V (T, p ) which was the average of character¬ 
istics with respect to jpr Q: 


V(T,p) 





( 6 ) 


Note that the dimension of the denominator in front of the integral and the 
dimension of integration limits cancel themselves out. 

From the topological properties of the characteristics’ set and on © as 
well as on © the pseudo-state equation for soft magnetic composites has been 
derived. This equation has enabled us to determine the optimal values of the 
technological parameters jb]. However, the described optimization relates only 
to the power losses and in design processes optimization of induction is also 
important. Therefore, the goal of this paper is to create a new algorithm for 
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optimization of both the power losses and the induction within the frame of a 
bicriteria problem. 

2. Experimental Data 

Specimens were produced by cold pressing under pressure of 500.. .900 MPa. 
The specimens made of Somaloy 500 powder were cured at a temperature of 
400.. .600°C for 30 minutes in air atmosphere. The specimens used in exper¬ 
iments were ring-shaped with a square cross-section. The specimens had the 
following dimensions: external diameter 55 mm, internal diameter 45 mm and 
thickness 5 mm. Total power loss density Ptot , expressed in watts per kilogram 
(W/ kg), was obtained from measurements of the AC hysteresis cycle according 
to IEC Standard 60404-6 using the system AMH-20K-HS produced by Labora- 
torio Elettrofisico Walker LDJ Scientific. Total power losses Ptot were measured 
at maximum flux density B m = 0.1... 1.3 T over a frequency range of 10 to 
5000 Hz. During measurements of the total power losses Ptot, the shape factor 
of the secondary voltage was equal to 1.111±1.5 %. Maximum measurement 
error of the total energy losses was equal to 3%. In order to optimize the mag¬ 
netic properties, the magnetic inductions B at fixed magnetic field H equal to 
1000 A/m were determined. These values were obtained from measurements of 
the DC magnetization curve according to IEC Standard 60404-4 using the same 
measuring system. 

3. Power Losses and Induction Pseudo-State Equations 

Optimization of the power losses was based on the topological properties of 
the characteristics. However, in the case of magnetic properties the situation 
is much simpler. For optimization of magnetic properties we selected induction 
Piooo f° r the fixed magnetic Held P=1000(A/m). We chose this value because 
the magnetic permeability of the soft magnetic composites reaches a maximum 
value around this magnetic field. We expected that the pseudo-state equation 
would properly describe induction at this point as a function of T and p. In 
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the previous paper [5| it was assumed and confirmed that the loss measure 
V obeys the scaling. Here, this assumption was extended to induction. In 
order to justify this assumption we referred to two phenomena: invariance of 
power losses (area of the hysteresis loop) with respect to scaling and invariance 
of the hysteresis loop with respect to scaling 7]. Therefore, for the bicriteria 
optimization problem, minimization of the power losses and maximization of 
the induction for a fixed magnetic field, we used the following pseudo-state 
equations of general form: 


where 


v(T, P )=^y . $(x), 

B 1000 (T,p)= yy ■ A(X'), 


X = 


x' = 


T_ 

Tc 


y Pc J 

T_ 

T' 

w 


(7) 

( 8 ) 

(9) 

( 10 ) 


where <!>(•) and A(-) were arbitrary functions to be determined. 7 , S , 7 ', S' 
and T c , p c , T' c , p' c are scaling exponents and scaling parameters respectively, and 
were to be determined. In the case of the power losses’ pseudo-state equation 
all calculations concerning modelling of $(•) and fitting of scaling exponents as 
well as model parameters were done in jfj] . The most important result was the 
derivation of an infinite set of solutions for the technological parameters which 
minimized the power losses: 

T_ 

7 ^ = 19 , 75 . ( 11 ) 

'Pc' 


4. Induction Pseudo-State Equation 

In this Section we derive a pseudo-state equation for induction H 1000 which 
will constitute a function of the two variables p and T. This function and the 
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Table 1: Somaloy 500. Measure of induction -Biooo vs. hardening temperature T and com¬ 
paction pressure p for magnetic field H —1000 (A/m). 


Temperature 

Pressure 

Induction 

(K) 

(MPa) 

m 

723,15 

800 

0,378 

773,15 

900 

0,496 

773,15 

700 

0,483 

673,15 

800 

0,335 

773,15 

600 

0,467 

823,15 

800 

0,546 

773,15 

500 

0,414 

741,15 

764 

0,425 

773,15 

750 

0,489 

773,15 

800 

0,504 

773,15 

650 

0,469 

773,15 

725 

0,467 

873,15 

800 

0,568 


power losses’ pseudo-state equation 0 will enable us to optimize induction and 
losses together. The optimization criteria are the following: find V = V m i n and 
Bwoo = -61000 max with respect to p and T. Deriving in Q the form for $(■) 
we reveal two phases of Somaloy 500: low losses and high losses. Therefore, in 
terms of the induction pseudo-state equation we have to take into account this 
phase separation. Measurement data of B 1000 vs. T and p are separated into 
these two phases in Table [1] 

The horizontal line between B woo = 0,414(T) and 61000 = 0,425(T) in¬ 
dicates the crossover between the low losses phase and the high losses phase. 
This transition is clearly visible in the jump of the V(T,p) function around 
the separation line js] . For each phase we assume an independent branch of the 
pseudo-state equation in the form of the Pade approximant. In order to simplify 
notations we introduce the following abbreviations: 


P 

TV — -, TV 

Pc 


/ 


P 


( 12 ) 


Expressing A(-) in @ by the Pade approximant we get the following form for 
the induction pseudo-equation of state: 


Biooo(T,p) 


/V Gp + Gi X' + G 2 X' 2 + G 3 X' 3 + G 4 jf ' 4 
1 + Lb X' + D 2 X ' 2 + D 3 X' 3 + D 4 X' 4 ’ 


where Go,..., G 4 , D±, ..., D 4 are parameters of the Pade approximant. 
parameters have to be determined from the data presented in Table [l] 


(13) 

All 

The 
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Table 2: Somaloy 500, low-loss phase. Values of the induction pseudo-state equation’s param¬ 
eters and the Pade approximant’s coefficients of CEL _ 


Y 

s' 

T' 

Pr. 

Go 

Gi 

g 2 

1,114 

0,499 

32,186 

25,849 

- 784,41 

764,05 

- 276,06 

g 3 

g 4 

Di 

d 2 

D 3 

Di 

- 

43,412 

- 2,4315 

3,6486 

2,9005 

2,8373 

3,975 

- 


Table 3: Somaloy 500, high-losses phase. Values of the induction pseudo-state equation’s 
parameters and the Pade approximant’s coefficients of CEL _ 


Y 

s' 

T c 

Pc 

Go 

Gi 

g 2 

1,1146 

0,4992 

32,19 

25,83 

- 808,91 

747,44 

- 266,95 

g 3 

Gi 

Di 

d 2 

D 3 

Di 

- 

42,852 

- 2,5142 

0,5846 

- 3,001 

- 9,968 

7,1432 

- 


corresponding pseudo-state equation for the power losses has been derived in 

d 


V(T,p) = ^ 


G 0 +G 1 X + G 2 X 2 + G 3 X 3 + G a X 4 
1 + D x X + D 2 X 2 + D 3 X 3 + D a X 4 ' 


(14) 


5. Estimation of Parameters for Induction Pseudo-State Equation 


The above-mentioned crossover between low-loss and high-loss phases is 
observed as A sudden change of V between two points: [773,15; 500, 0] and 
[742,15; 764,0] Tabled! However, this effect is not seen in the induction mag¬ 
nitude. Therefore in order to have a compact description of the power losses 
and the induction we take that into account and we divide the data of Table [T| 
into two subsets corresponding to the two respective phases. Minimizations of 
X 2 for both phases have been performed with Microsoft Excel 2010, where 

2 


N 


X' 


= s iooo(r/,7r') - 7if 


■Gq + Gx X 'j + G 2 X'X + G 3 X'\ + G A X'l 
1 + D 1 X'i + D 2 X ’ 2 + D 3 X'l + £>4 X'\ 


(15) 


where N = 7 and N = 6 for the low-losses and high-losses phases, respectively. 
Table [2] and Table [3] present estimated values of the model parameters for the 
low-loss and high-loss phases, respectively. 
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Figure 2: Scaled -Biooo vs. scaled temperature in the low loss and high loss phases. 


6. Optimization of Induction and Power Losses 

In the optimization of the power loss problem jjj we have applied low loss 
phase solutions and high loss phase solutions have not been considered. How¬ 
ever, it is not clear whether this simplification excludes important solutions for 
the induction. The binary relations are invariant with respect to scaling 3,2]. 
This enables us to present all scaled characteristics in the one picture Figj2] 
and draw the following conclusion. All considered pressure characteristics of 
the high losses phase are covered by the set of the low losses phase characteris¬ 
tics. Therefore for further investigations we limit our searching to the low losses 
phase. To this end we draw part of the phase diagram of Somaloy 500 corre¬ 
sponding to the low losses phase Fig. [3] and we deliver values of Gi,D i7 p c , 7 
which are displayed in Table [4j 
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Table 4: Somaloy 500, low-losses phase. Values of the V pseudo-state equation’s parameters 
and the Pade approximant’s coefficients of {HD 0. This table is an amended version of Table 
3 in 0 _ 


7 

5 

T c 

Pc 

Go 

G i 

g 2 

1,2812 

0,1715 

21,622 

37,729 

370315315 

-47752251 

1734952 

g 3 

Gi 

-Di 

d 2 

d 3 

Di 

- 

-1,3764 

-678,26 

170,80 

6243,8 

386,96 

-28,699 

- 


Table 5: Somaloy 500, low-losses phase. Optimum solutions in technological and in physical 
spaces. 


p( MPa) 

n° c) 

V(W/kgT“ p ) 

Biooo(T) 

389 

370 

14,1 

0,300 

492 

407 

20,0 

0,356 

584 

440 

27,3 

0,400 

683 

478 

40,0 

0,449 

733 

499 

50,0 

0,479 

764 

515 

58,5 

0,500 

800 

532 

70,0 

0,525 

838 

549 

82,9 

0,550 

906 

570 

101 

0,580 

979 

584 

116 

0,600 
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Figure 4: Technological optimum curve presenting dependence of optimum temperature vs 
optimum pressure. 
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Figure 5: Physical optimum curve. Induction B 1000 vs.losses measure V. 
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7. Details of Bicriteria Problem 


All calculations in this section have to satisfy the following conditions: 18,4 < 
X < 22, 9, which results from limitation of the presented calculations to the Low 
Losses Phase presented in Fig. [3] The considered bicriteria problem is formu¬ 
lated by the initial value of V = V\, feasible set of (p,T) and the two criteria : 
V{p, T) = V m i n whereas B WO o(p , T) = B W oomax ■ Since increase of .Biooo causes 
increase of V these conditions are in contradiction. Therefore looked for solving 
criterion should lead to self- consistency between V rnln and Biooo max- Such con¬ 
sistency will be achieved as a fixed point of the following recurrence procedure. 
Let maximiztion and minimization procedures be reprezented by the operators 
Omax and O m i n , respectively. Let Bi 000 (T,p) and V(T,p) be functions defined 
by m and m respectively. Then the one step of independent optimizations 
of B iooo (T,p) and V(T,p) can be writen in the following form: 

O max B1000 = B1000, max f OT T = T\, P = Pi (16) 

OminV = Vmin for T = T 2 , P=P2- (17) 

The obtained result consits of the two points (T\, p f), (T^i ^ 2 ) and the corre¬ 
sponding values of magnitudes to be optimized -BmooCTbPi) and V{T 2 ,P 2 )- 
Therefore, any further optimization is not possible and the bicriteria problem 
is not solved. In order to mesh B 1000 ,max and V m i n we introduce constrain 
V (T, p) = Vo which protects (TTfil) and (fTTl) against collapse, where Vo is an ini¬ 
tial value of loss. Then (fTTH) get the following form which coupled Brnoo, max and 


Vmin as well as leaved some space for further optimization: 

V(T,p) = V 0 , O max Bmoo(T,p) = B1000(Ti,pi). (18) 

Having Biooo{Ti,pi) we protect (fTTl) against collapse: 

Biooo(T,p) = B W oo{Ti,pi) 

5 Omin V = V(T 2 ,p 2 ). (19) 

Therefore, after n steps we obtain: 

V(T,p) = V(T 2n ,p 2n ), O max Bmoo(T,p) = B1000(T 2n+ i,p 2n+ i), (20) 

Biooo(T,p) = Biooo(T 2 n +l,P 2 n+l), OminV = V (T 2 n+2 , p 2 „ +2 ) ■ (21) 
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m-m generate the two converging series: Ti,T 2 , ■ ■ ■ ,T 2n +2 arid p \, p 2 , ■ ■ ■ ,p 2 n+ 2 - 


lim T k = T* (22) 

k—> oo 

lim p k =p*. (23) 

k—t oo 

Substituting T* and p* to (fl3l) and (UTl) we derive the meshed values of V and 
B 1000 : 

V* = V(T*,p*), B * 1000 = B WO o(T*,p*). (24) 

The found solutions are not unique. Selecting set of initial values for V 0 we de¬ 
rive the set of final solutions. Optimization has been done by SOLVER routine 
of EXCEL2010 program. Obtained output is presented in TABLE [5] Fig0]and 
FigE] present these results in technological and in physical spaces, respectively. 
The obtained results represented by markers are fixed points of the proposed 
procedure. There is one to one correspondence between these points in physical 
and technological spaces (TABLE [5]). In order to select an unique solution one 
must provide an additional criterion resulting from a relation between impor¬ 
tance of losses and induction. For instance, assuming the deepest minimum 
for the scaled measure of losses V7 r -7 we apply condition given by m • Inter¬ 
section of two curves presented in Fig [6] leads to the following single solution 
p = 382(MPa) and T = 363(°C). In the physical space this point corresponds 
to V = 13, 64(Wkg -1 T -/3 ) and Riooo = 0,29(T). 

At the end we pay some attention to the power losses measure V. This is an 
auxiliary magnitude which help us to derive values of designing technological 
parameters due to the following features: 

• V is pseudo-thermodynamic average with respect to magnitude created 
with the peak of induction and the frequency of electromagnetic field wave. 
Therefore, this includes information about both independent variables. 

• V depends on the technological parameters. 

• Physical dimension of V is unknown yet due to a dummy exponent (3. 
However, the value of V is well determined together with values of p and T 
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Figure 6: Reduction of the feasible set of solutions for the technological parameters {p, T} to 
the single point (p = 382MPa, T = 368° C). 


which enables us to compose SMC specimen and to perform measurements 
of its characteristics. Finally, applying 0,0 we are able to calculate /3 
and to determine the physical dimension of the current V. 

8. Conclusions 

We have presented method for the bicriteria optimization of the chosen phys¬ 
ical properties of Soft Magnetic Composites. By this way we have solved the 
problem mentioned in jfj] which concerns optimization of losses and induction. 
Achievement of the fixed point is interpreted as revelation of an equilibrium be¬ 
tween the both assumed criteria.The crucial roles in the presented method play 
scaling and the notion of pseudo-state equation. The created system is as good 
as the experimental data which have been used for the estimations of model 
parameters. Therefore, presented here the first version will be improved by 
forthcoming new experimental data. The presented example in this paper is a 
minimum nontrivial case of Multiphysics problem and shows that this approach 
suits for designing Magnetic Composites. Therefore, the presented algorithm 
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is going to be extended for more than two physical features of the composing 
material. For instance, the designing of magnetic composites requires also op¬ 
timization of mechanical properties, since the susceptibility of such materials 
to cracking in service is of fundamental concern Q. We address the derived 
algorithm to designers of SMCs. 
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